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Summary—The concept of additivity of individual resistances to mass transfer between phases 


is confirmed by experimental measurements of 


transfer coefficients for vaporization of water 


into air, and the absorption from air of carbon dioxide, ammonia, and acetone by water. 


Résumé—Les expériences de l’auteur vérifient la régle d’additivité des résistances individuelles 


aux échanges de matiére entre phases. 


Les coefficients d’échange ont été mesurés pour la vaporisation de l'eau dans lair et pour 
l'absorption par l'eau de gaz carbonique, gaz ammoniac ou acétone mélangés a lair. 


The industrially important of gas 
absorption and liquid extraction involve mass 
transfer between two phases, and equipment 
design procedures commonly apply the well- 
known “ two-film ” or * two-resistance ” theory 
{1}, [2]. Although this practice has been general 
for many years there would appear to be inade- 
quate experimental support for the concept of 
additive resistances, and it may be asked if, in 
fact, over-all coefficients can be calculated from 
values of the individual or * film ” coefficients. 
The usual theory assumes that there is no 
diffusional resistance of importance at the phase 
boundary, and the relationship between the con- 
centrations of diffusing solute in the two phases 


processes 


at the interface is assumed to be the same as in 
static equilibrium measurements with bulk phases. 
Several investigators have suspected the existence 
of a diffusional resistance caused by the concen- 
tration of high-molecular weight solutes at the 
phase boundary, but such an effect has not been 
demonstrated in mass transfer tests where the 
interfacial area was fixed and definite. Even 
though there is no diffusional barrier at the inter- 
face, the relationship between the concentrations 
in the two phases at the interface may be expected 
to be different under conditions of dynamic 





equilibrium with diffusion than under the condi- 
tions of static equilibrium for which the equili- 
This effect, which 
might be expected to indicate an apparent inter- 
face resistance varying with diffusion rate, has 
been shown by ScuraGe [3] to become important 


brium data are obtained. 


only at quite high rates of mass transfer, ordinarily 
attainable only at reduced pressures. 

The additivity concept has been tested recently 
in liquid systems by Gorpon and Suerwoop [4], 
who also review the question of additivity of 
resistances in transfer between phases. A large 
number of solutes were transferred between water 
and iso-butanol in a stirred vessel, under condi- 
tions of known interfacial area. The measured 
values of the over-all coeflicients were found to 
be well correlated using single values of the 
individual coefficients and the additivity concept. 
These results represent an excellent confirmation 
of the two-resistance theory, since the solutes 
used covered a 7,600-fold variation of the distri- 
bution coefficient. However, the values of the 
individual coefficients used in this correlation 
deviated 13 and 44 per cent, respectively, from 
the individual coefficients as measured directly 
by the extraction of iso-butanol into water and 
water into iso-butanol, with no third solute 
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present. This discrepancy was incompletely 
explained. 

The present work is a study of the addivity 
principle in gas-liquid systems. Again a stirred 
vessel was used, with a constant known area of 
contact between phases. The tests involved the 
vaporization of water into air, the absorption of 
carbon dioxide, ammonia, and acetone from air by 
water, and the absorption of ammonia from air 
by aqueous solutions of sulphuric acid. The test 
of the additivity concept involved the comparison 
of the observed coefficients for the 
absorption of ammonia and acetone with values 
predicted from the individual coeflicients obtained 
from the experiments on water vaporization and 


over-all 


carbon dioxide absorption. 


EXPERIMENTAL PROCEDURE 


The thermostated cylindrical glass absorption 
vessel shown in Fig. 1 was used in all tests. It 
was 14-2cm id. and 30-5cem high.  Air-gas 
mixtures were supplied to the upper gas space at 
constant condition and flow rate, but the liquid 
remained in place throughout each test. <A 
down-flow propeller agitated the liquid; the gas 
agitators consisted of the two flat-blade stirrers 
shown. Gas entering through a 14 mm o.d. glass 
tube was directed horizontally and tangentially 
to the blade tip of the lower gas stirrer, which 
barely cleared the end of this inlet tube. Calli- 
brated glass thermometers were used in both gas 
and liquid. Laboratory compressed air was used 
as the carrier gas for the various solutes, the rate 
of flow being constant at 1,512 + 201. per hour 
in all tests. The air was cleaned by being passed 
through a large bottle of glass wool and dried over 
calcium chloride. It was then metered, humidified 
by passing through a _  constant-temperature 
bubbler and a glass wool trap, preheated (or 
cooled) to the desired inlet test temperature, and 
passed to the absorption vessel. Heated air was 
used in the water vaporization tests ; in all other 
cases the inlet gas temperature was 25-07 + 
0-18°C. The liquid temperature was 24-95 4 
0-15°C. (In the sulphuric acid tests the liquid 
temperature rose from 25°C to about 26°C in the 
course of an experiment.) Test conditions (ex- 
cept liquid concentrations) remained exceedingly 
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constant. The inlet and outlet gas streams were 
sampled by aspiration and both analyzed by 
appropriate techniques. Liquid samples with- 
drawn for analysis were replaced by equal volumes 
of distilled water. 

In the case of the water-vaporization tests, the 
inlet air was preheated to a controlled dry-bulb 
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Fig. 1. Elevation of reaction vessel. 1, Thermostating 
bath; 2, vessel cover; 3, absorption vessel; 4, gas inlet 
tube; 5, liquid phase thermometer; 6, brace rod; 7, liquid 
replacement tube; 8, liquid phase stirrer; 9, supporting 
collar; 10, stirrer shaft; 11, liquid sampling tube; 12, lower 
gas phase stirrer; 13, upper gas phase stirrer; 14, gas phase 
thermometer; 15, mercury seal; 16, gas outlet tube; 17, 
access hole. 


temperature (30 to 45°C) such that its wet-bulb 
temperature was 25°C, This was done in order 
that the temperature of the water surface would 
be essentially the same as the bulk liquid tem- 
perature, and permit the use of the latter in 
calculating the vapour pressure at the phase 
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boundary. In these tests the partial pressure of 
water in the inlet air was varied from 0.00745 to 
0-02542 atm. in seven tests ; the partial pressure 
of water in the outlet air 0.0155 to 0.0272 atm. in 
the same series. Rate of vaporization was 
calculated from the air flow rate and the increase 
in partial pressure of water in the air. Air 
humidities were obtained gravimetrically by 
aspirating samples of some 161. of gas through 
drying tubes containing ‘Anhydrone’ (alumina was 
found to adsorb carbon dioxide from the air). 
Each test lasted approximately three hours ; 
data for only the middle third of the period were 
used. Explorations of the temperatures at some 
sixty points in a horizontal plane 13 mm above 
the liquid surface showed minor variations and 
a mean within 0-5° of the gas outlet temperature. 
It was concluded that the gas space was well 
stirred, and the outlet humidity was used in 
calculating the driving force for mass transfer. 

In the tests with carbon dioxide the inlet air, 
held at 25°C, contained 2 to 10-3 mole per cent 
carbon dioxide and was humidified to 98 per cent 
relative humidity to eliminate cooling of the water 
surface by vaporization. Tests varied in length 
from two to five hours. Carbon dioxide was 
supplied to the air stream after humidification 
from a pressure cylinder. Aspirated inlet and 
outlet gas samples were analyzed gravimetrically, 
using adsorption trains of ‘ Ascarite ’ and * Anhy- 
drone’ in series. Only one inlet and one outlet 
gas sample were aspirated over the duration of 
each test. Because the rate of absorption was 
very low, the compositions of inlet and outlet 
gas streams were nearly identical. Rates of 
absorption were calculated from the liquid 
analyses, made by titrating 100 ml. samples taken 
at approximately 30 minute intervals. The carbon 
dioxide-water equilibrium data of MacINNes and 
Be.cuer [5] were used in calculating the over-all 
driving force from gas to bulk liquid. 

The technique employed with ammonia was 
very similar to that for carbon dioxide. Both gas 
and liquid were at 25°, and the gas was nearly 
saturated with water. Inlet and outlet gas 
streams were analyzed for ammonia by titration 
following absorption of aspirated samples in 
bubbler tubes containing standardized hydro- 





Gas concentration was varied from 
0-3 to 1-1 mole per cent ammonia. Liquid samples 
were withdrawn into standard acid and _ back- 
titrated with standard sodium hydroxide. Gas 
samples were taken over approximately thirty- 
minute intervals. Liquid samples were withdrawn 
Partial 
pressures of ammonia over the week aqueous 


chlorie acid. 


at the end of each of these periods. 


solutions encountered were taken from the corre- 
lations and data of Suerwoop [6] and Breiren- 
BACH [7]. 

The procedure using acetone was also similar, 
except that acetone was not supplied from a 
cylinder, but acetone 
bubbler in the Aspirated gas 
samples were drawn through bubblers containing 


from a_ thermostated 


air feed line. 
distilled water. These solutions, as well as liquid 
samples from the absorption vessel, were analyzed 
by means of a Beckman spectrophotometer, 
Model DU, using a wavelength of 2,640 A. This 
procedure, which was found to be simple and 
accurate in the range of 0-2 to 6 wt. per cent 
acetone, was checked against solutions of known 
concentrations. 

Two tests were run in which ammonia was 
absorbed by aqueous solutions of sulphuric acid. 
The technique was essentially the same as for 
absorption of ammonia by water, except that the 
inlet air was humidified only to the extent of 
95 to 99 per cent of the vapour pressure of water 
over the acid solution used. 


Discussion or REsuULTs 


The technique employed in studying the vapori- 
zation of water, whereby the wet-bulb temperature 


of the air over the liquid was adjusted to be the 
the 
elimination of heat transfer into or out of the 
water, 


same as water temperature, ensured the 
and made it possible to assume with 
confidence that the water surface temperature was 
the same as that of the bulk liquid. The mass 
transfer coefficients obtained in the vaporization 
of water into air, therefore, are individual values 
of k,, not K,. Carbon dioxide is sufficiently 
insoluble in water so that the gas side resistance 
calculated from the observed k,, for water vapour 
is less than 0-01 per cent of the over-all resistance. 
The observed coefficients for carbon dioxide, 
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therefore, are taken to be values of k, rather 
than K,. 

Ammonia and acetone were chosen for the 
remaining tests, since their solubilities were such 
that calculations based on the observed values 
of k, for water vapour and k, for carbon dioxide 
indicated that for these gases, in the equipment 
employed, both gas and liquid resistances should 
be appreciable fractions of the total. 

The experimental results are summarized in 
Table 1, which shows the observed values of Kg 
or K, and outlet gas concentration in each 
test. Observed coefficients are seen to vary but 
slightly from the average in each series, though 


Table 1. Summary of data* 





( Og Absorption 
Outlet Partial 
Pressure, atm. 

0-0614 
0-0842 
00-1032 
O14 
01-0508 
0-0370 
(0260 
OO1eT 


Water Vaporization 
Outlet Partial 
Pressure, atm. 

0-0232 

0-0198 

O-O155 

00-0166 

00-0272 

0-000 

O-251 


Kg 
0-212 
0-200 
0-192 
0-204 
0-104 
0-204 
O22 
Ave. 0-208 
Range of CO, concentra- 
0 to 0-0017 


Range of vapour pressures 
0-03076 to 0-03098 atm. tions in liquid 
g moles /1. 
Ammonia Absorption Acetone Absorption 
00-0150 

ol 

00-0005 

O-O0TS 

0-0059 

o-0025 


00-0045 
OO077 
0O-0077 
00-0104 
OO1lO 
00-0085 
00-0058 
0-0041 
0-0030 
0-0045 
0-0060 


Ave 


Ave. 1-69 
Range of NH, concentra- 
tions in liquid 0 to 0-23 g 
moles /1. 


Range of acetone concen- 
trations in liquid - 0 to 0-11 


g moles /1. 





* The complete data may be obtained by writing to the 
second author. Copies of the twenty pages of original data 
will be supplied for the cost of reproduction. 
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gas concentrations were varied appreciably. Liquid 
concentrations varied from zero to 17 x 10*¢g 
mole /l. for carbon dioxide, zero to 0-23 g mole /I. 
for ammonia, and zero to 0-11 gmole/l. for 
acetone. One acetone run was discarded, since 
the value of K, was 12 per cent below the average 
for the other runs. 

The test of the concept of additive resistances 
consists in comparing the observed values of K, 
for ammonia and acetone with the values calcu- 
lated from the observed coeflicients for water 
vapour and carbon dioxide. In making this 
comparison it is assumed that k, is proportional 
to Dj, and k, proportional to D7, where D,, and 
D, are the molecular diffusivities in air and 
water, respectively. Fortunately, D, for 
ammonia is very nearly the same as for water 
vapour, and D, for ammonia is very close to 
the value for carbon dioxide. Consequently, the 
comparison of observed and calculated values of 
K, for this solute depends very little on the 
assumed effects of diffusivities. 

Many investigators in the field have employed 
the power function to represent the effect of 
diffuisvity, using values of n and m in the range 
of 0-15 to 0-67. Both exponents will be assumed 
to be 0.5, since this is not only representative of 
much of the experimental data, but is the value 
suggested by the penetration theories of Hicare 
[8] and of Danckwerts [9]. The addition of 
resistances then gives : 


1 l 
, b, Os 
= (dD, 1.0 


H (ke, 


M,, 


(kveo.| cp, + oO, 


| 


Using values of D, and D, for water, carbon 
dioxide and ammonia from Perry [10], D,. for 
acetone from the GiLLILaND [11] equation, and 
D, for acetone from WiiKr [12], the following 
ratios were obtained : 


Acetone 
0-367 
0-712 


Ammonia 
0.9021 
1-17 


Dy, (Dd, Mo 
D,. (Dico, 
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The comparison of observed and calculated 
coefficients is summarized in Table 2. Using the 
diffusivity ratios given above, the calculated 


Table 2. Summary of Results 





Calculated 
Kg Ky 


Experimental 
Ke Ay, 
Vaporization of water 0-203 
Absorption of CO, 
Absorption of NH, by water 
Absorption of acetone by water 
Absorption of NH, by acid 
(normality = 3) 


3-02 

1-o 

1-95 
O-1s82 





value of K, for ammonia is 1-76, which is to be 
compared with the observed value of 1-69. The 
calculated and observed coeflicients for acetone 
are 1-89 and 1-93, respectively. The liquid side 
resistance is 45 per cent of the total in the case 
of ammonia, and 62 per cent for acetone. Thus 
the observed within 
3-5 per cent with the values calculated from the 
water and carbon dioxide data. The relative 
unimportance of the diffusivity correction in the 
case of ammonia is indicated by the fact that 
the predicted K, is 1-73 when no diffusivity is 
made. The data lend excellent support to the 
concept of additivity of individual phase resist- 


over-all coeflicients agree 


ances, and suggest that over-all coeflicients may 
be calculated with assurance from the individual 
coefficients. The concept of a substantial resistance 
at the interface appears to be incompatible with 
the data, except for the unlikely possibility that 
such a resistance, expressed in liquid-side con- 
centration units, is the same for ammonia, carbon 
dioxide, and acetone, except as modified by the 
diffusivity correction. 

The ammonia data give no insight into the 
effect of diffusivities since the ratios are so near 
unity. For acetone, however, the effect of the 
lower values of D, and D,, should be significant. 
Assuming » and m equal, it is found that the 
observed and calculated values of K, agree if 
these exponents are taken as 0-47. This is very 
close to the value 0-50 so commonly used in the 
literature. 


ABSORPTION IN ACID 


The absorption of ammonia by acid solutions 
should presumably follow the Harta theory [2], 
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with Kg increasing as the acid strength is 
increased until K, becomes equal to k, and 
then remaining constant as acid normality is 
increased above some critical value. The data 
of Fe.urncer [18] show an anomaly in this 
connection, however, since K,a for the absorp- 
tion of ammonia from air by sulphuric acid in a 
packed tower was found to increase steadily 
with acid strength up to 10N, and reach values 
more than twice as great as the apparent ka as 
calculated from the K,,a for ammonia absorption in 
water by allowing for the relatively small liquid- 
side resistance. In order to explore this point, 
two tests were run in which aqueous sulphuric 
acid was used for the absorption of ammonia. 
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4 
ACID NORMALITY 


Absorption of ammonia from air by water and 
aqueous solutions of sulphuric acid. 


Fig. 2. 


The results are shown in Fig. 2, with K, 
plotted vs. acid normality. Although the data 
are incomplete, they fit well qualitatively with 
the Hatta theory. The over-all coefficient Kg 
increases with acid normality, reaching an 
apparent asymptote of about 0-182, which is not 
far from the value of 0-192 calculated from the 
value of 0-203 for water, with allowance for 
diffusivity based on the square root relation. The 
results do not support the anomaly suggested by 
the Fellinger packed tower data, which was 
possibly due to the equivalent of a variation of 
a with acid normality. 
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NOTATION 


= gas-side coefficient, g moles /(hr) (em®) (atm). 
overall-coefficient, g moles /(hr) (em?) (atm). 
liquid-side coefficient, em /hr. 
over-all coefficient, em/hr. 
molecular diffusivity in liquid, em? /hr. 
molecular diffusivity in gas, em®/hr. 
Henry's law constant, atm/(g moles /1). 
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Summary 


may be altered fundamentally by different methods of feeding the slurry to the basket. 


July 1952) 


It has been demonstrated that the mechanism of cake formation in a hydroextractor 


¢xamples 


are given of cakes growing in a direction parallel to the basket wall rather than normal to the 
wall, leading to the formation of wedge-shaped cakes. 


Résumé 
de coin, aussi bien pour l’essorage. 


Etude théorique et expérimentale de écoulement au travers d'un gateau en forme 
L’auteur compare l'equation précise de rélaxation a des 


relations simplifiées (écoulement du liquide uniquement dan la direction du champ de forces). 
Discussion d’une equation basée sur l'emploi d'une épaisseur uniforme, équivalente a l’épaisseur 


moyenne du coin. 


Attempts to formulate equations representing 
the formation of hydroextractor cakes will depend 
on the mechanism controlled by the operating 
conditions and characteristics of the packing 
material. A cake may be formed with varying 
thickness by spraying slurry on to a cake fed 
with filtrate at less than the critical wetting rate 
[1], the solid remaining where it first hits the 
cake. Cakes of fine materials formed while sub- 
merged beneath a slurry layer may become either 
a wedge shape or a cake of nearly uniform thick- 
ness, depending on the characteristics of the 
material and the feed arrangement. Ideally, a 
slurry fed uniformly to the face of a hydroextrac- 
tor basket builds a cake of uniformly increasing 
thickness at all heights on the basket wall. With 
a feed to the bottom of the basket various 
materials are found to build preferentially from 
the base of the cake and unless remoulded the 
final cake is wedge-shaped [5]. 

The final wedge shape is not necessarily pro- 
duced by proportional growth over the whole 
surface of the cake. The cake may be formed 
predominantly from the bottom, the upward 
growth giving a thinner cake near the top. Thus 
the following limiting mechanisms could be for- 
mulated mathematically and compared with 
experimental facts : 


(a) a cylindrical cake grows evenly at all posi- 
tions, 





(b) the cake grows upwards between two speci- 
fied radii. 


Case wirh RapiaL Growrn or CAKE 

As one limiting mechanism consider the formation 
of a cake submerged beneath a slurry layer, i.e. 
r,<7,. Assume that the slurry is fed to a 
horizontal basket at a rate to maintain r, con- 
stant, and that the cake grows evenly at all 
positions with equal increments of thickness. If 
the formation is begun on a cylindrical face, 
either of cloth or of cake material, the cake face 
remains vertical and cylindrical. 








i= 
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Hypothetical types of growth for derivation of 
equations (5) and (10). 


Fig. 1. 


When an aqueous slurry of concentration s is 
fed to a hydroextractor basket spinning at n and 
the cake grows towards the axis of basket, keeping 
a smooth cylindrical surface (Fig. la), the cake 
surface radius r, is related to the volume V of 
the slurry added from the start of formation : 


* Present address : Institutute of Technology, Alexandria University, Alexandria, Egypt. 
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sV 


7,3 — ——— 
a-X-op 


. (1) 


where C, 
The rate of feed of slurry, g.. 


* 


drainage of liquid through the cake, q,. 


& 7 X Tp: 

and the rate of 
are 
related by : 


It has been shown [1], [2] that for a liquid layer 


from r, to r, on a cake running full of liquid 


L 
from r, to ry, : 


— nn? KX (7? — 1,7) 


TA 
; 

ug log, -“ 
r 


c 


This equation was appropriate to the cakes 
tested as the sedimentation rate in the nominal 
slurry layer from r, x 
layer was effectively liquid only. 


to r, was so high that the 


When a layer of thickness (r, — r,) and per- 
meability A is formed on an initial layer of 
thickness (7, — 7.) and permeability A,, the flow 


rate of liquid through the banded cake is 


tr? in? XC 7. > 
a ug (r,° : ry) ; | 
(Sa) 


» }-1 
. log, | 


0 'o , 


when forming cakes on highly permeable cloth 
K, may be fifty times A and with (rv, — r,) 

(rg —7,) the second term in the bracket is 2°, 
of the first term and thus the cloth term soon 
becomes of minor importance when building a 
thick cake. It is questionable, however, whether 
this type of equation is valid once the low per 
meability layer is appreciable, as the highly 
permeable cloth will probably cease to run full 
of liquid, and atmospheric pressure is trans- 
mitted to the radius r, [4]. In view of the short 
time required to reach this stage, the equation 
(3) has been used, which accords with the assump- 
tion of a large value for K, in equation (3a). Thus 
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CX 
q, — 
r 8 
log, —. ug « (1 - 
r o 


c 


r,?) 


2K (,2 
where C 4n° . A (To 


CX 
iT 
‘ log, ( . ro. . 

(7? Cc, 3 ) 
Oatt Otolv =v 
assuming C and C, to be constant, 


When integrated from V 
at ft i, 
equation (4) becomes : 


(1 


+ 2 log, 7) V 


log, (r,? —€ 2CAT 


, 2r,? 
iF) c log, Ty 


1 
Case with Urpwarp Growin or Cake 
If the cake is assumed to grow in a direction 
parallel to the axis of rotation and confined 
between radii r, and r, (Fig. 1b), the relation 
, and V is given by 
sV 


Op 


between ry? 


nx (r,* — r,*) (6) 


‘ . 9 v 
Cc, Ve. where C, 


map (1,\" — 73") 


Assuming that this cake is being formed on an 
existing cylindrical layer of inner radius r, and 
the same permeability K, 


tr® n? K (7.2 — r,*) 


me-(r—2) 


whence 


R, V+ R,, where R, 


R, — log, 
" 


From integration with V = 0 at t = ¢, assuming 


R, and R, to be constant, 


it 
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(9) 


(10) 


The significance of the equation is limited to the 
range 0 <a < X. 

Fig. 2 is a comparison between the theoretical 
equations (5) and (10) for the set of experimental 





~— 


eo 


w 





——-— equation. 5S. 


a 


—— equation.lO. 











o 
> 
; 
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34 
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is 2025 
TIME, t minutes. 


10 3% 


Fig. 2. Comparison between equations (5) and (10). 


constants given in Table 1. The time required 
for forming a certain cake according to equation 
(10) is much less than that to be predicted from 
equation (5), especially when the building starts 
on a thin layer. The calculations required for 
the use of equation (5) are extensive and require 
seven figure accuracy. 

Curve A is from equation (5) based on a 
negligibly small initial resistance. Curves B, C, 
D are from equation (10) representing different 
initial cake layers from r, to r,. The meet of 
curves B and A represents the volume of slurry 
that formed the initial cylindrical layer, which 
could have been formed by the mechanism of 
equation (5), the formation following equation 
(10) thereafter. 





Hydroextraction VIII : Cake formation 


EXPERIMENTAL RESULTS 

A twelve litre glass vessel was equipped with a 
mechanical agitator, a feeding tube to the floor 
of the basket and an inclined graduated tube 
from the side for measuring the level of the 
slurry in the vessel. Dilute starch slurries of 
selected concentration were used for the tests, as 
chalk cakes were found to crack occasionally 
under the test conditions. 

The hydroextractor basket was fitted with two 
layers of cloth support of high permeability and 
run at a specified speed with water being fed to 
it to maintain the liquid inner radius constant 
at a selected value of r,. This was done using 
the probe arrangement as an indicator for r, 
[2], [3]. The slurry was then fed to the basket, 
starting at the moment when the water feed was 
stopped, and the rate continuously regulated to 
maintain r, constant. The level on the calibrated 
side tube was recorded every 15 or 30 seconds 
according to the rate of feeding. 

Concentrations of 0-10 to 0-20 g air-dried 
starch /ce slurry were used. Cake formation was 
studied by feeding slurry directly on to the cloth 
support or on to an existing cylindrical cake of 
smooth surface, and known dimensions and 
permeability. Tests were done at different speeds 
of rotation and different radii r, and r,. Constant 
dimensions and cake characteristics appropriate 
to representative tests are quoted in Table 1. 

The surfaces of the resultant cakes and of 
intermediate forms were traced using the probe. 
To give a clear visual indication of the way in 
which the cake grows, in some tests the addition 
of slurry was interrupted at intervals to feed a 
small amount of a solution of starch and iodine. 
After completing the formation, the cake was 
run dry for some minutes and then vertical 
sections cut from it and withdrawn from the 
basket. Fig. 3 shows sections of representative 
cakes along the plane of the axis of rotation, the 
contours of which were measured by the probe 
[2] for comparison with stained sections. It 
demonstrates that the cake formation starts with 
a wedge at the bottom of the basket. This 
increases in thickness up to the inner radius of 
the slurry r,. At this stage, the wedge may be 
of more than an inch thick at the bottom with 
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3. Contours representing the growth of hydro- 


extractor cakes. 


zero thickness at the top of the basket. Thereafter 
the main growth of the cake is in a direction 
parallel to the axis of rotation, with hardly any 
This 


mechanism suggests that the experimental V’, ¢ 


increase in cake thickness at the bottom. 


data will not agree with the function given from 
equation (5). The data in Fig. 3 show a similar 
different 
different positions of the liquid surface, and for 
formation on cloth (Figs. 3a, b) or on an existing 
cake of constant thickness (Fig. 3c). 

Curve A in Fig. 4 is predicted from equation (5), 
curve B from equation (10) for the set of experi- 
mental results represented by the data, using the 


formation at speeds of revolution, 


data in Table 1 representing the known geometry, 


and characteristics of the cake material, The set 


A, B is for a cake formed on an existing layer of 


starch, whereas the set of curves A’, B’ is for a 
cake formed directly on the cloth. If the cloth 
is assumed to have negligible resistance to flow, 
equation (10) becomes inapplicable. Hence 
curve B’ has been made to represent a cake which 
grew from an initially uniform thickness of one 
millimetre. This is expected to lead to a slightly 
higher time than the exact value that could be 
obtained from equation (10) alone taking the 
resistance of the cloth into consideration, if the 
cloth was assumed to remain full of liquid. 

In both cases, the failure of equation (5) is 
clear, whereas equation (10) gives theoretical 
nearly with the 
results. A variety of cases is represented in Fig. 5 
where the agreement between the 


lines coincident practical 


practical 
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Fig. 4. 
curves. 


Comparison of practical results with theoretical 
The arrow indicates the volume of slurry equiva- 
lent to the cake volume to fill the basket. 


results and predicted curves from equation (10) 
The data all the 


curves are presented in Table 1. 


is confirmed. relevant for 


In the above calculations the overall mean 
permeability, including the effect of cloth-cake 
interface resistance [3], was considered constant, 
and taken as the permeability of a cake formed 
from half the The 
permeability actually increases from outer to 
inner layers of the cake [3] as well as from bottom 
to top of the basket. 

Correction for the former type of variation 


material used in the test. 


will not produce a significant change in the 
deviation of equation (5) from the experimental 
data. In the case of equation (10), the latter 
variation in permeability would decrease the 
curvature of the predicted curve and make it 
more nearly coincident with the experimental 
data. However, this variation is of minor impor- 
tance when testing starch as the porosity was 
nearly the 
through the hydroextractor cakes. 


constant along vertical direction 
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Fig. 5. Agreement between practical results and equation 


(10). The arrow indicates the volume of slurry equivalent 
to the cake volume to fill the basket. 


CONCLUSIONS 


Even with submerged cakes, formation may not 
be assumed to occur by growth all over the cake 
surface, when a dilute slurry is fed to the base 
of the basket. By assessing face contours during 
stages of the formation, an approximate equation 
may be derived from the hydroextraction equa- 
tions, as in the cases described which conform 
closely to the upward growth mechanism, The 
appropriate equation depends on the materials 
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Table 1. 








op = 0-75 g/ce 
wn = 0-0114 poise 
to = 108 cm 


X = 11-1 em 


Test conditions and cake characteristics 


ao = 1-59 g/cc 













Ss 
g solid n’ 
| 

Curve Figure cc slurry T.p.m. 


1080 
1080 


1080 


1080 

















used and the operating characteristics. 


A 5 0-15 900 3-0 8-0 9-1 
B 5 0-10 1200 3-2 77 9-62 
Cc 5 O-15 1900 3-1 8-2 10-0 
D 5 0-20 2100 30 6 |)=— 677 10°5 
NOTATION 
C, Cy, Cg = constant groups 

g = gravitational acceleration cm/sec.” 

K, Ky= permeability g/sec.? 


n = speed of revolution 

q = flow rate, q, for liquid, q, for 
slurry 

R,, Rg = constant groups 

r = radius, r, for cake inner surface, 
ry, for liquid surface, rg, for cake 
outer surface, r,, for inner surface 
of basket wall 

s = slurry concentration; (g air- 
dried material/ce slurry) 

t = time 

)’ = volume of slurry 

X = depth of basket 

x = height of cake 

p = Viscosity of liquid 

ao = density ; o for air-dried solid ; 
op for packing density 
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Summary 


This paper gives binary equilibrium data for the miscible systems benzene-acetic 


acid, ethyl acetate-acetic acid, acetic acid-water, and for the partially immiscible system ethyl 


acetate-water. 


Different values of the end activity coefficients are obtained for ethyl acetate 


and water, depending on the method of calculation. Furthermore, it is shown that in the immis- 
cible region the extrapolated van Laar relationship does not hold. 


Résumé 
acétique ; 
incomplete : 
acétate d'éthyle 


acetate d'éthyvle — acide acétique ; 


acétate d'éthyle — eau. 


Données d'équilibre sur les systemes binaires enti¢rement miscibles : 
acide acétique 
Suivant la méthode de calcul, on obtient pour le systéme 
eau des valeurs différentes pour les coefficients d'activité limites. 


benzene — acide 
eau ; et sur le systéme a miscibilité 


En outre, 


dans la région dhétérogénéité, la relation de Van Laar ne peut étre extrapolée. 


INTRODUCTION 


The purification of acetic acid by straight distilla- 
tion is difficult due to the low relative volatility 
of 


concentrations. 


acetic acid-water mixtures at high water 
The separation of these two 
components has been simplified by employing 
solvent extraction followed by azeotropic dis- 
tillation. 
feeding into the top of the column a third com- 


the 


The azeotropic distillation consists of 


ponent which greatly increases relative 


volatility of water to acetic acid and at the same 


time is only partially miscible with water. 
Ethyl acetate has frequently been used com- 


mercially as the entrainer in this azecotropic 


separation but no vapour-liquid equilibrium data | 


for the ethyl acetate-acetic acid-water system 
could be found in the literature. It is essential 
for the accurate design of an azeotropic column 
that the ternary vapour-liquid equilibrium data 
should be available or be capable of prediction 
In this the 


prediction of the ternary equilibrium data from 


from the binary data. instance, 
the binary data is complicated duc to the associa- 
tion of acetic acid in the vapour state in both the 
binary and ternary systems. For this reason the 
binary and ternary vapour-liquid equilibria have 
been determined experimentally. The purpose of 
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this paper is to present the vapour-liquid equilibria 
of the binary systems investigated. 


Mernop or ANALYsIs 
The results on the benzene-acetic acid system 
were obtained using density determinations at 
25°C and by titration with barium hydroxide 
solution for acetic acid as a check. The density 
values obtained are given in Table 8 


Table 8 
Densities of Benzene- Acetic Acid Mixtures 





Mole %, Benzene Density at 25°C 
10442 
101193 
1-00614 
0-98873 
0-95132 
093338 
0-90226 
O-88522 


100-0 0-872 





The ethyl acetate-acetic acid samples were 
analyzed using similar analysis methods. The 
densities obtained for mixtures of ethyl acetate 
and acetic acid at 25°C are given in Table 9. 
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Table 9 
Densities of Acetic Acid-Ethyl Acetate Mixtures 





Mole %, 
Ethyl Acetate 


Density 


a 25°C 


1-0442 
1-01356 
0-98372 
0-96527 
0-93563 
0-8945 





The acetic acid-water samples were analyzed by 
determinations of the acetic acids by titration 
with barium hydroxide solution when less than 
60%, acetic acid by weight was present. If more 
than this amount was present the water content 
was determined by titration with Fischer reagent. 
The two methods were shown to agree over the 
range where both methods could be used. 

The ethyl acetate-water samples were analyzed 
for water content using Fischer reagent. 


Purity of REAGENTS 
Acetic acid, ethyl acetate, benzene, water, were 
purified and analyzed as described in Part 1 [19] 
and Part 2 [23]. 


EquiLiprium DETERMINATIONS 


The apparatus described by Exits [24] was used 
in most cases, using essentially the same proce- 
dure. In certain cases where indicated a modified 
GILLEsPIe [24] still was used. The temperature 
of the vapours leaving the still were maintained 
1-2°C above that of the boiling liquid. Samples 
were drawn off into ice cooled bottles after 
hours operation. 


Vapour Pressure Dara 


The vapour pressures of benzene were calculated 
from the equation : 
8-11401 


logy, P (mm Hg) 2046-74 


] 
7T( A) 


given by the National Bureau of Standards [25}. 
The vapour pressures of acetic acid were cal- 
culated from the equation 
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1206-35 


6-89745 aetna 
220(273 + t) 


log,, P (mm Hg) 
given by Brown and Ewan [21]. 

The vapour pressure data for ethyl acetate 
taken International Critical Tables 


were from 


, ] 
[26] and correlated by plotting log P ys. 1 


Activity COEFFICIENTS 
In this investigation activity coefficients have 
been calculated using the equation 


(1) 


where y, = activity coeflicient of component 1 
total pressure 

vapour pressure of component 1 at 
the temperature of the determination 
mole fraction of component 1 in the 
vapour 

mole fraction of component 1 in the 
liquid 


BENZENE-ACETIC AcID 


The results given in Table 10a were calculated 
using the molecular weight of monomeric acetic 


Table 10a 
Benzene- Acetic Acid 





TC 


111-2 
92-3 
89-3 
BOD 63-0 

54-2 
45-0 

B47 


25-0 





Table 10b were calculated 
using the molecular weights of acetic acid vapour 
at the temperatures given. The molecular weight 
data by Rirren and Simons [27] were used. The 
results are expressed graphically on Fig. 6. It 
will be noticed that the results based on Table 10b 
are in much better agreement with the van Laar 
equations than the data of Table 10a. 


acid. The results in 
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Fig. 6. 


Table 10b 
Benzene- Acetic Acid 





M.W. of 
Xp Acetic 
Acid 


22-9 
68-8 
74:5 
79-4 
82-1 
85-7 
(88-5 
840 91-0 


101-0 
103-0 
103-5 
104-0 
104-1 
104-4 
104-6 
104°8 


95-06 771 
67-4 31:3 
58:3 25-5 
48:55 20-6 
0-6 17-9 
3190 143 
234 115 
160 90 





The evaluation of the activity coefficient of 
acetic acid in solution with another component is 
made difficult by association of acetic acid in the 
vapour state whether as a pure component or in 
a mixture. If acetic acid is present in a vapour 
mixture, and the second component does not alter 
the extent of association of the acetic acid com- 
pared with that of the acid in the pure state at 
that temperature, the associated acetic acid may 
be treated as a pure component (i.e. as distinct 
from a mixture of polymeric forms). The calcu- 
lation of the liquid and vapour compositions is 
then made using the molecular weight of the 
acetic acid vapour at that temperature. If the 
second component in the mixture influences the 
state of association of the acetic acid in the 
vapour mixtures, the above procedure is no 
longer valid. The investigation of Raman spectra 
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Extraction of acetic acid from water: Part 3 
has shown that benzene does not influence the 
extent of association of acetic acid in liquid 
mixtures and it seems reasonable to assume that 
it will have no effect in the vapour state. It was 
thought, therefore, that if the molecular weight 
of the acetic acid in the pure vapour state at the 
relevant temperature be used to calculate the 
molar fractions in liquid and vapour mixtures of 
benzene and acetic acid, the activity coeflicient 
curves will obey one of the equations given by 
Caritson and Co_nurn [13]. This procedure has 
been suggested previously by Burier [28] and 
by Herincron [29). 

In the comparison (Fig. 6) of experimental and 
van Laar activity coeflicient plots, temperature 
and gas law deviation effects have been neglected. 


Eruy. Acetrate-Acetic Acip 
The results obtained for this system are given in 
Tables lla and 11b. together with the molecular 


Table lla 


Ethyl Acetate- Acetic Acid 





x. 


97-8 “67 
85-4 ‘92 
73-2 1-05 
58-9 1-06 
42-2 1-07 
32:7 1-07 


12-8 
43-9 
58-0 
98-5 
S18 
86-3 





Table 11b 


Ethyl Acetate- Acetic Acid 





M.W. of 
Acetic 
Acid 


X4 Yep 


105 
103-5 
102-8 
102 
101-5 
101-1 


1-29 98-7 


79-5 
88:°3 
93-6 
99-7 
105°1 
108-1 


7:74 
32-6 
44-6 
59-5 
72:8 


79-0 


90-7 
82: 

70-9 
55-3 


45-0 


9-4 
17-7 
29-1 
44:5 


55-0 





weights used for acetic acid in calculating the 
mole fractions. The log activity coefficients given 
in Table 11b are plotted in Fig. 7. It can be seen 
that the deviation from ideality in this system 


is very small and consequently small errors in the 
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Fig. 7. 


analysis of the vapour or liquid compositions 
markedly influence the log activity coeflicients. 
Difficulty was experienced in analyzing mixtures 
by titration for acid when the acid was present 
in high concentrations. The results were checked 
by determining the density of the mixtures. Fig. 7 
suggests that probably the best correlation of the 
experimental points is given by the two Margules 
curves. It seems possible that a certain amount 
of inter-association may take place giving curves 
which differ slightly from those in Fig. 7. 


Acetic Acip-WaTEeRr 
The acetic acid-water system has been investi- 
gated quite extensively [21], [30], [31] with a 
variety of equilibrium stills. 


Table 12a 


Acetic- Acid Water 





1 1005) 827 91-7 663 93-4 1-43 1-00 1-64 
2 101-3 19-7 803 13-6 86-4 1-20 1-03 6-12 
3.1022 2902 7OR8 201 799 +116 1-04 9-1 
$ 103-2 38-7 61-3 278 727 116 105 109 
5 106-1 61-0 39-0 48:5 51-5 117 106 12-5 
6 1070 66-0 340 53:3 46-7 1:16 1-07! 12-65 
7 108-6 73-9 26-1 61-6 384 1-14 109 12:3 
8 100-7 77-9 221 66:1 33-9 1:12 1:09 118 
9 111-8 85:6 144 760 240 1-10 1-11 9-6 
10 1128 91-5 8&5 843 15-7 107 1-16 7-2 
11 1160 965 355930 705 1-04 3-5 
12 104-7 52-0 48-0 39-5 60-5 1-17/106 12-5 
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little consistency between the results obtained, 
it was decided to reinvestigate the system using 
the equilibrium still mentioned above. 


Table 12b 
Acetic Acid Water 





102 


1 1005 S02 95-0 40 960 1:42 995 

2 101-3 12-6 87-4 85 91-5 1:17 1-00 101-8 

3. 1022196 80-4 12-9 87-1 I-11 1-00 101-8 

$# 103-2 27-1 72-9 185 81:5 1:10 1-00 101-7 

5 106-1 51-5 485 381 61-9 1:09 1-03 101-3 

6 1070 53:5 465 404 596 1:08 1-00 101-2 

7 108-6 626 37-4 489 S11 106 1-01 101 

8 109-7 67-6 324 536 464 105 1-02 101 

9 1118 780 220 654 346 104 1-04 100-9 
10 113-83 865 13:5 764 23-6 103 1-09 100-7 





The results obtained are listed in Tables 12¢ 
and 12b, 
data with that obtained by previous investiga- 
tions a plot of y,-7, for water was made. These 


For the purpose of comparison of the 


values are given in Table 12a and plotted in 


Fig. 8. The results obtained in this investigation 
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As there has been | 
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agree well with those of ConneLL and Monronna 
[30] and are in reasonable agreement with those 
of Brown and Ewatp [21]. The results of Gar- 
win and Happap have a slightly higher relative 
volatility in the middle portion of the concentra- 
tion range but agree well with the present data 
at the extremities. 

The log activity coefficients for this system are 
plotted in Fig. 9. Correction for the molecular 
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Fig. 9. 


weight of acetic acid in the vapour still does not 
give a reasonable correlation. This is probably 
due to the change in the degree of association of 
acetic acid when diluted with water as pointed 
out by Traynarp [31]. No correlation of this 
system was attempted in the absence of quanti- 
tative data on the change in association of acetic 
acid vapour in the presence of water vapour. 


Eruyt Acerate-Water System 

The activity coefficients of the two components 
in the ethyl acetate-water binary system are 
considerably higher than those for the ethyl 
acetate-acid and water acetic acid systems. In 
the prediction of ternary equilibrium data in the 
ternary ethyl acetate-acetic acid-water system 
from those of the binary systems the activity 
coefficients of the ethyl acetate-water system are 
therefore the controlling ones. Unfortunately, 
this binary system is the most difficult to investi- 
gate experimentally. 

As may be seen from Fig. 10, the ethyl! acetate- 
water system presents the particular difficulty of 
being immiscible in both the liquid and vapour 
condensate phases over the greater part of the 
concentration range. Any liquid concentration 
falling within the range AB will give a vapour 
of composition represented by D. This means 
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that, once this azeotrope composition has been 
determined, only the vapour-liquid equilibria of 
liquids with compositions between A and E and 
B and C have to be investigated. The analysis at 
each end of the composition range is thus liable 
to be inaccurate. The composition at A for the 
ethyl acetate-water system is only between 1 and 
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2 mole %. A further problem concerning the 
equilibrium determinations is that over a con- 
siderable portion of the above concentration 
ranges, the vapour condensate will be two-phase 
due to the composition being within the mis- 
cibility boundary FA BG. 

Table 13 lists the vapour liquid equilibrium re- 
sults obtained for this system and these terminal 
results are plotted on Fig. 11. The composition 
of the azeotrope was found to be 69-7 mole per 
cent ethyl acetate which agrees well with other 
results in the literature. 


Table 13 
Ethyl Acetate-Water System 





3 
~ 


Yp  Xw YE 


Yw 


1-02 
1-02 
1-02 
1-02 
1-05 


91-4 | 286) 
04-7 1-69 
95-4 1-41 
92-7 | 2-59 
95-3 | 1-57 


eo oe 
a oS @ 


Ge 
-_ =~ 
—— 


~~ 7 -) +) +7 





In all, three methods were used for the pre- 
diction of end values (i.e., log y, at X, = 0, 
log yp at X,, = 0 etc.) of the activity coefficients 
of this system. These were as follows : 


(a) by experimental determination of the vapour- 
liquid equilibria when the liquid and vapour 
condensate are both single phases ; 

(b) from the azeotrope composition and the com- 
position of the two liquid phases in equili- 
brium. When the boiling liquid mixture is 
partially immiscible, the activities of these 
phases will be equal and consequently each 
liquid phase will be in equilibrium with the 
azeotropic vapour composition ; 

by using the boiling points of the miscible 

liquid mixtures between 0-1 mole °, and 

90-100 mole %, of ethyl acetate and extra- 

polating these activity coeflicients to find the 

end values. 


Griswo._p et al. [32] have investigated the 
vapour-liquid equilibria of the ethyl acetate-water 
system and obtained an end value of 8-0 for the 
activity coefficient curve of water. In the eventual 
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correlation, however, GriswoLpD used the solubi- 
lity equations given by CoLBuRN and ScHEON- 
BORN [33] and end values of 6-5 and 76 were 
obtained for water and ethyl acetate respectively. 

CaRLsOoN and CoLBuRN have stated that the 
solubility equations are known to give results 
which vary considerably from those obtained 
from vapour liquid equilibrium determinations. 
This arises from the assumption that in the 
original derivation of the solubility equation the 
van Laar or Margules equations are valid. 

The activity coefficient end values for water 
obtained in the present investigation are given 
as follows : 

8-3 
8:3 


(a) Vapour liquid equilibria data 
(b) 
(c) 


Cc 


From azeotrope determinations 


From boiling point measurements 11-0 


The results obtained from boiling point data 
have been found to be higher than those obtained 
from equilibrium determinations. This would be 
expected since the method of calculation is based 
on an approximation. Thus the figure of 8-3 
would appear to be best of the end values for the 
activity coefficient of water in this mixture. 

The end values obtained for ethyl acetate are 
as follows : 


48-0 
65-4 
50-0 


(a) From vapour liquid equilibria 
(b) From azeotropic data 
(c) From boiling point data 


The vapour liquid equilibrium results for ethy! 
acetate in water were determined on a still used 
by SELVEPATNAM [834] in which both the immis- 
cible liquid and vapour condensate mixtures were 
stirred. 

The above activity coefficient end values have 
been mainly determined the activity 
coeflicients of miscible mixtures which exist up to 
the extremities of the two-phase boundaries and 


from 


of these values the first is preferred. 
In considering the immiscible 
problems arise : 


region two 


(a) Can the binary van Laar curves in the 
miscible region be extrapolated into the 
miscible region and expressed in terms of the 
overall composition of the liquid phases ? 
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(b) Can the binary van Laar curves be extra- 
polated in the increased miscible region pro- 
duced by the presence of the third compo- 
nent ? 


HERRINGTON [29] has stated that within the two 
phase regions it is possible to use activity coefli- 
cients based on the overall composition of the 
liquid. Assuming this is correct then the activity 
coeflicients in the immiscible region are defined 
as follows : 
log K, (2) 
! a 


log 


| 7 
og P 


Y 
4 


log y2 = log "* 


le mu 
Por, ” a 


(3) 


K, and Ky are constants since in the immiscible 
region the temperature is constant and a hetero- 
geneous azeotrope is formed giving constant values 
of y, and y, for binary composition x, and x,. Thus 
if the temperature and azeotrope composition is 
known, A, and A, can be evaluated. 

Equations (2) and (3) are of a different form 
from cither the van Laar or Margules equations 
and consequently when plotted on Fig. 11 for 
ethyl acetate-water the curves for equations (2) 
and (3) do not coincide with the extrapolated 
van Laar curves ; therefore it would appear that 





the binary van Laar curves in the miscible 
region cannot be extrapolated into the immiscible 
region. 

From Fig. 11, it is interesting to note that the 
activity coefficients for the boundary compositions 
evaluated from equations (2) and (3) are lower 
than those evaluated from the extrapolated van 
Laar curves in the miscible region. Whilst this 
difference within experimental error, 
Ropinson and GiuuiLanp [35] also refer to the 
deviation from the van Laar equations of experi- 


may be 


mental data of systems approaching immiscibility 
and suggest that this is due to entropy effects 
occurring on mixing. Thus, until further experi- 
mental work is available, question (b) cannot be 
satisfactorily answered; this information is, 
however, essential for the accurate prediction of 
vapour liquid equilibria in the immiscible ternary 


region. 


NOTATION 


P = vapour pressure 
T = temperature K 
{ = temperature “C 
a = total pressure 
¥, = mole fraction of component 1 in the vapour 
2, = mole fraction of component 1 in the liquid 
y = activity coefficient 
K,, Ka = constants 
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Liquid-liquid extraction in a ‘‘ rotating disc contactor” 
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( Received 19 February 1954) 


Summary—Using the system water-acetic acid-methyl isobutyl ketone, the extraction per- 
formance of a “ rotating dise contactor ~ (R.D.C.) of small diameter was investigated for various 
values of the rotor speed, total throughput and solvent to feed ratio. By comparing the separating 
efficiency with the fractional volume of the dispersed phase under the same circumstances, it 
was found that under certain conditions the efficiency decreases, although the hold-up of the 
dispersed phase increases. This effect is ascribed to back-mixing in the continuous phase due 
to entrainment by the dispersed phase. 


The effectiveness of the R.D.C. is compared with that of similar extraction columns. 


Résumé—En employant le systéme de leau-acide acétique-méthylisobuty'cétone extraction 
& contre-courant dans un “ rotating dise contactor” (R.D.C.) a été examinée pour plusieurs 
valeurs de la vitesse rotatoire, du débit total des liquides, et du rapport entre les débits des deux 
liquides. Une comparaison entre l'efficacité de séparation et le volume fractionnel de la phase 
dispersé sous les mémes conditions montre qu’en certains cas l’efficacité diminue pendant que la 
teneur en phase dispersée augmente. Ceci est attribué au mélange longitudinale dans la phase 
continue dai a lentrainement par la phase dispersée. 


L’efficacité de 'R.D.C, est comparée a celle des colonnes d’extraction semblables. 


INTRODUCTION are related to the interfacial area between the 


During the past few years several new types of | 'W° phases. 


extraction apparatus have been described which In this work no attempt was made to investigate 
show a greater efliciency and a better operational the effect of the dimensions or the constructional 
flexibility than the more conventional types details of an R.D.C. on the separating efficiency. 
(sieve-tray column, packed column, spray tower). Neither were the conditions investigated under 
In particular we refer to extraction columns which which Hlooding occurs. The extraction experiments 
have a rotating element : the Scuemper extractor | WET restricted to or system only, the well-known 
[1], [2]. the * rotating core column ” [3], [4], the system water-acetic acid-methy! isobutyl ketone 


(MIBK). 


* multi-stage mixer column” [5] and the “ ro- 
tating disc contactor” [6]. In these extraction 
columns, the intimacy of the contact between the EXPERIMENTAL Drralis 
two liquid phases depends primarily on the rotor | (a) The liquid system 

speed. By adjustment of this new independent | ‘There were several reasons for choosing the system 
variable the most favourable conditions for a| water-acetic acid-MIBK : ‘ 


given extraction process can be achieved. = , ' , , 
5 I 1. The acid content can easily be determined by 


titration with NaOH to the phenol phthalein 
endpoint. 


To determine these conditions the influence of 
the main operating variables, i.e. rotor speed, 


total throughput and solvent to feed ratio, on the Sgt ; 
extraction performance of a rotating dise con- | ?: rhe distribution ratio of acetic acid between 


tactor (R.D.C.) was investigated. At the same water and MIBK is near unity, and depends 
time the relative hold-up of the dispersed phase but little on the concentration of the acid or 
was determined under various conditions. A the temperature. 

comparison of the column efficiency and the The system shows little tendency to form stable 
hold-up was made since both of these quantities emulsions except after prolonged use. 
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. Since this system has been used by many other 
investigators a comparison of the performance 
of different types of equipment is possible. 


Equilibrium measurements had previously been 
carried out in this laboratory. The results were in 
close agreement with those of ScuerBe. and 
Karr [2], which have been used for further 
calculations (see Table 1). For the densities of 
both phases, which depend on solute concen- 
tration, the data published by OLpsuur and 
Rusuton [5] were used. 

Table 1. 


Equilibrium distribution of acetic acid between water and 
MIBK at 20°C. 





Acid content of MIBK-phase 
(2, kg /m*) 

Acid content of water-phase 
(y, kg/m*) 


0 100 200 30-0 400 


0 200 38:5 55-8 71-4 





In all experiments to be described acetic acid 
was extracted from MIBK with acid-free water. 
Before each run the water- and MIBK-phases 
were saturated with respect to each other. After 
the extraction, the resulting extract phase con- 
taining water, some acid and a small amount of 
MIBK discarded. The raflinate phase 
(MIBK) was recirculated for further use. 

The acid content of the MIBK feed varied 
between 3 and 5 per cent by weight. At these low 


was 


concentrations both the equilibrium line and the 
operating line are sufficiently straight to justify 
the use of a rapid numerical calculation of the 
efliciency. Furthermore, the physical character- 
istics do not change appreciably during the 
extraction process. Only the interfacial tension 
cannot be sufficiently controlled; even at low 
acid contents it depends markedly on the con- 
centration. Moreover, its value is affected by the 


inevitable small amounts of impurities. 


(b) The extraction column 

Fig. 1 shows the construction and the dimensions 
of the “ rotating dise contactor’ It con- 
sists of a vertical cylindrical tube fitted with a 
number of horizontal flat rings, spaced at equal 
distances. The rotor, which is mounted coaxially, 
is provided with a number of equally spaced flat 


, 


used. 





discs. 
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Sectional view of R.D.C. 

















All dimensions in mm 














MIBK +ACID | 


9..x5 


"Water + ACID 
$o.¥e 
Fig. 1. The extraction column (R.D.C.). 

As has been pointed out by Reman [6] in his 
original paper on the R.D.C., the turning of the 
rotor generates in the liquid in each compartment 
two torcidal vortices, which rotate in the hori- 
zontal plane. In this way a rather uniform degree 
of turbulence is produced, the intensity of which 
is governed by the rotor speed for given dimen- 
sions and liquid system. It is believed that by 
this arrangement an emulsion of relatively uniform 
drop size is produced since spots with excessive 
turbulence or very high rates of shear are absent. 

The top and bottom parts of the column contain 
the bearings of the rotor and the feed and exit 
lines for both liquid phases. Top and bottom 
were provided with an after-settler as indicated 
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in Fig. 1, with return lines for the heavy and light 
phase respectively. In most runs these settlers 
appeared to be superfluous, since nearly all of 
the settling took place in the free space in the 
top and bottom of the column. Separate experi- 
ments proved that the amounts of acid extracted 
in the settlers were negligible. 

In view of the corrosive nature of the system 
used, the column was constructed of stainless 
steel. The bearings were of a self-lubricating 
material, Philite. Neoprene rubber was used as 
a packing material. 
this had to be renewed several times during the 
year. 


Because of some swelling, 


(c) Auviliary equipment 
The feed and 
solvent were fed from high reservoirs by gravity 


This was of conventional nature. 


and the incoming flow rates were controlled by 
means of stainless steel needle valves and measured 
with calibrated rotameters. All runs mentioned 
here were performed with the heavy (water) phase 
dispersed in the light (MIBK) phase. By adjusting 
the bottom discharge stream (either manually or 
automatically with an on-off controller responding 
to changes in dielectric constant) the interface 
between the two liquid phases in the lower 
settler was kept at a constant level. 

For reasons mentioned above, all parts were 
made of glass or stainless steel. For junctions 
between glass and steel lines natural rubber 
tubing, which had to be replaced from time to 


time, was used, 


PROCEDURE AND CALCULATIONS 


(a) Extraction 

Before starting a run the entire column was filled 
with MIBK, which contained the acetic acid to 
be extracted (concentration 2,), and the MIBK 
flow rate was fixed at a desired value. After that 
the acid-free water was fed into the top of the 
column at a constant desired flow rate. A steady 
operation was obtained after 15 to 30 minutes. 
Attainment of the steady state was checked by 
rapid calculation of the material balance of the 
acetic acid from the observed flow rates, x,, and 
the acid contents of the MIBK- and water-phase 
leaving the extractor (7, and y, respectively). In 





general this material balance checked within 2 
to 3%,. All experiments were carried out at room 
temperature, i.e. between 17 and 22°C. 

The efficiency of the extraction operation was 
calculated in terms of the product of Kp, the 
overall mass transfer coefficient with reference to 
the dispersed (water) phase, and the total area of 
interfacial contact A. Since the equilibrium and 
the operating lines were approximately straight, 
K,A_ could the 


equation : 


be calculated from simple 


Pyy, — KyA (Ay) » 


where 


(y,* —y%) — y,* 


( Ay \im. 


In [(yo* — yp)/y*) ” 


The values of y,* and y,* (in equilibrium with 2, 
and a, respectively) were taken from a graph of 
the equilibrium relationship between 2 and y. 

The use of equation (1) is justified provided 
there is no back-mixing in either phase. The 
values of K,A are mean values for the whole 
column under the conditions prevailing at the 
time of the measurement. 

The column efficiency has also been evaluated 
in terms of the “height equivalent to one 
(H.E.T.S.) according to the 
method of VaRTERESSIAN and Fenske [7]. For 
these calculations the effective height of the 
column was taken to be 0-47 m. 


theoretical stage ”’ 


(b) Holdup 


For the determination of the fractional volume of 
dispersed phase (¢,) several methods can be 
used, For instance, at a certain moment all feed 
and exit lines can be closed simultaneously, after 
which the fractional volume the 
column can be measured. This method does not 
give the true average ¢,, of the extraction section 
since the liquid in the top and bottom parts is 
included in such measurements. 

In principle ¢, can be found by determining 
the average density of the liquid in the extraction 
section by amanometric method. Results obtained 
in this way were not sufficiently accurate because 
of the short length of the column and because of 
wetting difficulties ; moreover, there were small 


of water in 
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additional pressure effects due to the swirling 
motion of the liquid in the column. 

The method consisted of 
suddenly draining a part of the contents of the 
extraction section and determining the fractional 
water content after this sample had settled. To 
this end the column had been fitted with a short 
side-tube near the lower end of the extraction 
section (see Fig. 1). By suddenly opening a 
wide-bored stopcock in this tube during normal 
operation of the extractor, an amount of 150 to 
250 mil. (about 1.3 of the total volume) was 
drained from the extraction section within a few 
Smaller samples reduced the precision 
if samples larger than 


finally adopted 


seconds. 

of the oy observed : 

250 mil. were taken, a part of the continuous 

phase below the extraction section was entraind.e 
These 

always carried out in duplicate. In 75°, of the 

cases the two determinations checked within 5°. 


relative hold-up measurements were 





Discussion or ReEsu.ts 

Two series of experimental runs were performed : 
one in which the solvent to feed ratio was kept 
constant (®,/®,. = 1), and the total throughput 
was varied (12 <(®, + ®,) < 60l/hr), and a 
second in which ®, /®,. was varied between 0-25 
and 4 at a constant total throughput (®, + @,.) 
of 30! hr. In both series four different rotor 
speeds were used (350, 600, 1200 and 1800 r.p.m.), 
The efficiencies and ¢,-values are reported in 
Tables 2 and 3. 


(a) Efficiency 

Fig. 2 shows that K,,A is approximately propor- 
tional to the total throughput if the solvent to feed 
ratio is kept constant. The value of KA at a con- 
stant throughput, however, passes through a maxi- 
mum when the value of ®, /®,. is varied (Fig. 3). 
This maximum is located at a value of ®,, @, 
between 1 and 2 when @,, + ®, = 301 hr. 


Table 2. 


Extraction efficiency and fractional volume of dispersed phase in R.D.C. (1.D. = 41 em) Op)/ Oo = 1. 





Rotor Water MIBK 
speed 


rpm. 


350 
350 
350 
350 
350 
600 
600 
600 
600 
ooo 
Hoo 
600 
1200 
1200 
1200 
1200 
1200 
1800 
1800 
1800 
1800 
1800 


MIBK 


Water 10? 10° x 
% KpA 
kg om" Mig ‘sec. 


H.E.T-.S. 


10 
1-8 
2-6 
3-9 
5-2 
12 
2:3 
3-3 
42 
43 
5-2 
60 
18 
2-8 
47 
6-4 
8-9 
2-2 
42 
78 
6-9 
8-2 


170 
21-2 
20-5 
19-0 
22:1 
14-0 
25-6 
25-3 
141 
26-2 
21:8 
21-4 
19-5 
27-3 
27-09 


0-54 
0-46 
0-46 
0-48 
0-39 
0-36 
0-33 
0-30 
0-27 
0-26 
0-23 
0-18 
0-21 
0-19 
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Pp r 





Extraction efficiency and fractional volume of dispersed phase in R.D.C, (1.D. = 4-1 cm). 





, 






= B01 hr. 











No. Rotor MIBK MIBK Water 107 » 10° 
of speed ®1)/ Pc 2p x > bn KpA 
run r.p.m. kg /m® kg/m® kg /m® m® /sec 
23 350 0-25 38-1 30-4 30-3 31 1-1 
(3) 350 1-0 32-5 13-6 20-5 71 2-6 
24 350 20 38-6 7:48 15-7 8-6 2-9 
25 350 40 38-6 3-54 8-80 9-7 2-3 
26 600 0-25 37-4 29-0 36-0 34 1-4 
27 600 0-50 26-3 14-9 22:8 48 2-3 
(8) 600 10 36-0 11-9 25-5 6-2 3-3 
28 600 20 26-3 5-67 11-1 738 2-8 
“u 600 3-0 29-0 420 8-82 9-1 24 
bl 600 +0 29-0 2-80 7-06 o-4 24 
31 1200 0-25 27:8 18-4 36-4 +4 25 
32 1200 0-50 27-9 12-4 30-0 +S 3-7 
(15) 1200 10 35-2 8-05 27-9 8-6 4-7 
33 1200 20 B52 3-06 16-8 10-4 4-7 
oe 1200 20 45-2 371 21-0 oS +5 
35 1200 40 32-0 1-04 785 11-9 34 
36 1800 0-25 22-0 12°5 36°5 6-7 7 
37 1800 0-25 32-0 17-8 53-3 76 a+ 
38 1800 0-28 29-3 12-5 45-6 9-2 6-5 
39 1800 0-5 29-6 8-86 41-1 11-1 71 
w 1800 0-5 23-2 7-32 31-3 9-3 6-3 
(20) 1800 10 29-6 2-97 27-0 16-2 738 
(21) 1800 10 23-2 2-81 20-9 13-7 oo 
(22) 1800 10 39-1 3-65 35-1 16-7 8-2 
41 | 1800 2-0 29-3 1-00 14-7 19-2 O-+ 








HE.T-.S. 


m 


0-88 
0-61 

0-46 
O-45 
0-70 
0-60 
0-45 
0-48 
0-46 
O45 
0-45 
0-36 
0-33 
Ost 
0-32 
0-33 
0-23 
0-20 
0-20 
O-19 
O22 
0-18 
0-21 

0-19 
0-21 
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Fig. 2. K,A vs total throughput. 
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350 r.p.m. 
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1200 
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Fig. 3. A, vs. solvent to feed ratio. 


An increase of the rotor speed markedly | smaller. As a result of this, of course, the 
improves the efficiency of the column because the | maximum throughput of the column is limited 


average drop size of the dispersed phase becomes | to a lower value. In these experiments, for 


20 
O/b-=1 


Oo 350 r.pm 
= a * 
© 1200 

& 1600 








” 


” 




















40 


Fig. 4. Volume fraction of dispersed phase vs. total throughput. 
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instance, the maximum total flow rate (®,,+@,.) 
for ®,/®, = 1 was about 35 |/hr at 1800 r.p.m. 
and 501/hr at 1200 r.p.m. 

As can be seen from Tables 2 and 3, the 
H.E.T.S. values of the column are influenced 
primarily by the rotor speed. Under favourable 
conditions, i.e. not too small a liquid load and 
®,/®,. of the order of unity, the H.E.T.S. 
appear to be about 20, 30, 45 and 60 cm for 1800, 
1200, 600 and 350 r.p.m. respectively. 


(b) Volume fraction of dispersed phase 
The hold-up measurements reported in Tables 1 





Figs 4 and 5 show that ¢, increases with 
increasing rotor speed and with increasing flow 
rate of the dispersed phase, as was to be ex- 


pected. 


(c) Relation between K,A and ¢), 


Figs. 2 and 4 and also the left-hand side of 
Figs. 8 and 5 show great similarity which suggests 
a parallel between KA and the fractional volume 
of the dispersed phase. Therefore in Fig. 6 these 
two quantities have been plotted for various 
values of the rotor speed. A rather good propor- 
tionality between K,A and ¢p is observed as 





20 


| 
6,+0,.= 30 I/hr 


350 r.p.m 
600 





rs 
8 
s 


1200 
1600 














lj 
iP 














Fig. 5. 


and 2 were preceded by a preliminary investiga- 
tion concerning the influence of flow rates and 
rotor speed on ¢,. Throughout these preliminary 
runs the binary system water-MIBK was used in 
order to avoid variations in physical properties, 
in particular surface tension. Consequently these 
results showed much less scattering than the 
values reported for actual extraction conditions. 
The latter ¢, values were always somewhat 
higher than those obtained with the binary system 
under the same conditions. 
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—_ 6/9, [\/hr] 


Volume fraction of dispersed phase vs. solvent to feed ratio. 


long as ¢p does not exceed the value of 0-15 and 
®,/®, is not much greater than 1. 

As has been noted by Apper and Exern [8], 
SHERWOOD, Evans and Lonccor [9], Jounsox 
and Buss [10] and Grer and Hovcen [11], under 
certain conditions in spray towers there is a 
proportionality between K,A and ¢,. In these 
cases it is assumed that both these quantities are 
directly proportional to the number of drops of 
the dispersed phase per unit volume, with the 
mean drop size remaining constant. For spray 
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towers the latter will be approximately true, since 
for a given system the drop size is primarily 
determined by the nozzle construction. 

In the R.D.C. 


depends mainly on the flow conditions in the 


the drop size distribution 


extractor, which are determined primarily by 





10 
1800 rpm. 











oe 
| 





——> KA [10 ©m/sec] 


1200 r.p.m 


























15 20 
= 1009) 


Fig. 6. Relationship between A, 4 and fractional hold-up 

of dispersed phase. (— ——-—) = constant throughput 

and varying solvent to feed ratio ; ( )} = constant 
solvent to feed ratio and varying throughput. 


the rotor speed, so a proportionality between AA 
and ¢, can be expected for a constant speed of 
the rotor. 


Fig. 6 shows that this relationship between 
K,A and ¢, obtains if the hold-up is increased 
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by increasing the total throughput at a constant 
value of ®,/®,. of unity. 

If, on the other hand, ¢, is increased by 
increasing the solvent to feed ratio, ®, ®,, 
while the total throughput is kept constant, the 
The 
most probable explanation for this is that at 
high ®, ®, values the continuous phase is 


value of K,A passes through a maximum. 


entrained by the dispersed phase to a marked 
extent. This special kind of * back-mixing ” in 
the continuous phase impairs the efliciency of 
the extraction operation, This is reflected in the 
relatively low values of K,,A, as calculated by 
means of equations (1) and (2), since these 
equations are not valid when appreciable back- 
mixing occurs in one or both phases. Similar 
effects have been observed by Grer and Hovucen 
{il} OLDSHUE 
Rusuron [5}. 

In our experiments, at a value of ®,)/®,. > 1 


and been mentioned by and 


the back-mixing effect appears to become import- 
ant irrespective of the rotor speed and_ the 
hold-up. This conclusion should not be generalised 
however, since not only the solvent to feed ratio, 
but also the equilibrium relationship and the 
column efliciency itself determine for a certain 
type of column the conditions at which back 
mixing is significant. As the effect is of importance 
for the operation of most extraction columns with 
a high efliciency and a relatively small volume, 
further study is desirable, especially of the 


entrainment and subsequent mixing of the 


continuous phase. 


COMPARISON WITH OTHER EQUIPMEN' 


Reman [6] has pointed out that for a comparison 


of different types of extraction equipment both 
the efficiency (e.g. expressed in H.E.T.S.) and 
the flow of liquid handled should be taken into 
account. As an index of effectiveness of operation 
he proposed: J, = total throughput /volume of 
one theoretical stage. In Table 4 a number of 
values of this index have been collected for the 
acid-MIBK. As _ far 
possible, the most favourable examples were 
taken from the literature. Strictly speaking the 
comparison is not complete, if the value of 
®,,'®,, the choice of the dispersed phase, the 


system water-acetic as 
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Table 4. 


Comparison of effectiveness of different types of extraction equipment. 








Acid 
Type Dis- content Rotor Column H.E.T.S Effective- 
of Reference persed | Solvent ®,/%_ | in feed | speed, I.D., m ness 
column phase weight = r.p.m. m index I, 
% min™! 
Core [3] Table I No. 25 water water 0-83 6 1000 0-021 0-094 0-8 
Core [3] Table I No. water water 1-35 6 1000 0-021 0-15 0-7 
Core [3] Table I No. 53 water water 0-57 6 1075 0-012 0-029 1-9 
Scheibel [2] Table VII No. 75 | MIBK | MIBK) 1-0 15-20 400 O-287T | O-24 1-4 
Scheibel [2] Table VII No. 85 MIBK water 1-0 15-20 693 0-287 0-32 15 
Scheibel [2] Table VII No. 96 water MIBK 1-0 15-20 560 0-287 0-24 1-7 
Scheibel [2] Table VII No. 95 water water 1-0 15-20 400 0-287 0-32 1-5 
Scheibel [2] Table VII No. 100 water MIBK 1-0 15-20 400 0-287 0-23 2-1 
Scheibel {1} Table 12) No. 18 ? MIBK | abt. 1 20 1680 0-025 0-086 1-4 
Scheibel {1} Table 12. No. 19 ? MIBK | abt. 1 20 1680 0-025 0-056 2-2 
Scheibel [1] Table 12) No. 20 ? MIBK abt. 1 20 1680 0-025 0-061 2-0 
Mixer [5] Table 16 No, 12 MIBK MIBK 20 3 300- 0-153 0-14 2-5 
Mixer [5] Table 16 No. 4 MIBK MIBK 20 3 -450 0-153 0-099 2-0 
R.D.¢ [6] Fig. 8 water water 10 ? 1835 0-079 0-061 3:3 
R.D.¢ [6] Fig. 8 water | water 2-0 ? 1835 0-079 0-059 3-4 
R.D.C Authors not reported MIBK — water 20 3-5 1200 0-041 0-21 1-9 
R.D.¢ Authors | in this paper MIBK | water 10 3-5 1800 0-041 0-21 3-2 
R.D.4 Authors MIBK _ water 10 3-5 1200 0-041 0-25 1-6 
R.DA Authors Nos. 20-22 water water 10 3-5 1800 0-041 O19 | 2:1 
R.D.( Authors No. 41 water water 20 3-5 1800 0-041 0-21 1-9 
acid concentration and the direction of mass | J. M. pe Murxcx Keizer who participated in a 


transfer are different. Roughly it can be said, 
however, that there is not much difference between 
the “ multi-stage mixer column” [5] and the 
R.D.C., which in fact are very similar to each 
other. On the average the Scuerne. columns 
show somewhat lower figures; they have relatively 
small H.E.T.S. values, but their specific capacity 
is also lower than that of the two other types of 
extractors mentioned. 
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NOTATION 


A = total interfacial area in extractor 
I, = “ effectiveness index,” total throughput divided by 
volume per theoretical extraction stage 
Ky, = overall mass transfer coefficient with reference to 
the dispersed phase 
# = acid concentration in MIBK-phase 
y = acid concentration in water-phase 
¢p = volume fraction occupied by dispersed phase 
®- = continuous phase volume flow rate (here MIBK) 
®, = dispersed phase volume flow rate (here water) 
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conditions of viscous flow 
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Summary—tThe dynamic equilibrium of a suspension of uniform spherical particles settling 
in a fluid has been considered, and an expression has been obtained for the drag force on a con- 
stituent particle. This has been expressed in the form of a correction factor to be used in 
conjunction with Stokes’ Law, and has been compared with the correction factor obtained from 
experimental measurement. Two configurations of particles have been considered and good 
agreement has been obtained between the theoretical and experimental correction factors at 
volumetric concentrations greater than about 0-05, though the correction factor extrapolates 


to zero, corresponding to zero drag on the particle, at infinite dilution. The reasons for this 


di crepancy are discussed. 


Résumé—Etude de l’équilibre dynamique de la sédimentation dans un fluide d'une suspension 
de particules sphériques uniformes: établisement dune expression pour la friction. Cette 
derniére est donnée sous forme d'un terme correctif s‘ajoutant a la loi de Stokes; elle a été 
comparée au terme correctif déduit des résultats experimentaux. Les auteurs ont etudié deux 
configurations des particules et ont obtenu un accord satisfaisant entre les facteurs de correction 
theorique et experimental tant que la concentration volumétrique dépasse 0,05, bien que ce 
facteur puisse s’extrapoler jusqu’ a zero, correspondant ainsi 4 une friction nulle, pour une 
dilution infinie. Discussion des raisons de ce désaccord, 


In 1851 Strokes [1] obtained a theoretical 
expression for the drag foree (F) acting on a 
single spherical particle moving relative to a 
continuous fluid. 


Fk Ooo ps Vad (1) 


where « = viscosity of the fluid, 
d = diameter of the sphere, 
V’ = relative velocity. 


The fluid was assumed infinite in extent and | 


the relative velocity low so that the inertia of 
the fluid could be neglected. Equation (1) has 
been tested experimentally and found accurate 
at low velocities; for a Reynolds Group with 
respect to the particle of less than 0-2, the error 
does not exceed 4°). 

Since many industrial operations, such as 
sedimentation, elutriation, and _— tluidization, 
involve relative motion between a fluid and a 
system of suspended particles, it is important to 


be able to predict the forces acting on the par- 
ticles, and hence their settling rate in the fluid. A | 
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number of experimental studies have been made 
and empirical relations developed for the effect 
of concentration on the rate of sedimentation of 
suspensions of particles, but no satisfactory 
theoretical treatment of the problem has emerged. 


/ 


/ \ 


ee 





Fig. la. Hexagonal pattern, 
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In the present paper a number of simplifying 
made to enable a 


theoretical treatment to be applied to the problem 


assumptions have been 
of the sedimentation of a uniformly dispersed 
suspension of spheres under conditions where the 
resistance is due entirely to viscous drag. Because 
the 


forms a sharp interface, it 


and 
is assumed that all 


suspension settles at a constant rate 
the particles are settling at the same constant 


velocity and that a statistically stable arrange- 
ment is maintained. 

In order to simplify calculations, McNown [2] 
suggested that each particle in a suspension could 
be regarded as being situated in the centre of a 
Fig. la. The 
layers were taken 
particles 


Fig. Ib). 


hexagon of fluid as indicated in 
vertical distances between the 
as equal to the distance between the 
in a horizontal direction (Arrangement 1, 


ii 
WI 


Particle configuration (J) b/R 
10° 














7} 
| 








1126 °\/C at 
volumetric concentration. 


Fig. Ib. 


On the other hand, Hawks.ey [2] considered the 
particles to be arranged in adjacent horizontal 
layers so as to offer the minimum resistance to 
the motion of a fluid flowing through the system 
(Arrangement 2, Fig. 1c). In this second arrange- 
it is assumed that the hexagonal pattern 
in a horizontal plane is retained. The relation 
between the volumetric concentration of particles 


(C) and the ratio of the particle radius (6 = d/2) 


ment, 
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Particle configuration (11) b/R = 1-225 1/C at 
10°, volumetric concentration. 


















































Fig. lec. 


to the distance between the centres of the particles 


in a horizontal plane (2 R,) is then obtained as 


follows 


Arrangement 1 
Volume of each hexagonal prism = 4 4/3 R,* 


_$7h 


~ 44/3 BR! 
ae b \* 
24/3 (7) 


(2) 


Volumetric concentration = 


> _ 1193 4/6 
R, V 


Arrangement 2 


: 1/3 R,? b 


Volume of hexagonal prism 


Volumetric concentration 


j 226 ./T 
‘ R, 1-286 4/ ¢ 

In order to simplify the boundary conditions, it 
is assumed that each particle is surrounded by a 
cylinder of fluid with the same cross-sectional area 
as the hexagonal prism, so that the radius, R, of 
the cylinder is given by 


oR? = 24/3 KR? 


and 
R 
RK, 


1-05 
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For the first configuration 





b 1-183 , 
R 1-05 





VC = 1126 4/C 





and for the second : 









bh «1286 
R 1-05 V 











velocity. 





EQUATIONS OF FLOW 
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moving upwards at a velocity V, 


is equivalent to that of a suspension settling at a | 
velocity V, relative to a fixed horizontal plane. 
Let r be the radius of an elementary ring of the 
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The resistance force to the motion of a * 
sentative ” particle is dependent on the distances 
between the particles in the horizonatal layers. In 
the first configuration the particles are closer and 
therefore the velocity gradients are steeper and 
the drag is greater for a given mean relative 


For convenience we shall regard the particles as 
stationary and the free surface of the fluid as 





. € d | dV | » d j \ 
* Sep. 1 (dl) (r =— J dn Qa. 1, Ldn) (rP) dl 


| r) 
(4) 


d | 
(dl) Gn 


on On 







\ 


r Pe 





= (dn) dl 


(3) lan . (rP) dl 
AL 


fluid round a particle, dl be the thickness of the | i.e. 


ring in the direction of the streamlines, and dn 
the width along the equipressure lines as shown | 


in Fig. 2. 


Vis the velocity of the fluid along a streamline 
at point (r, h) and @ is the angle between the 


streamline and the vertical. 


Fig. 2. 
particle, 


EQuiLiprium or THE RING 


Since the ring of fluid is in equilibrium, the 
resultant shearing force balances the force due 


to the difference in fluid pressure, P. 
Thus 


av 
(2rrdl. . pe) dn 
i 


on A d 


Diagram of elementary ring of fluid round a 


‘ (27 rdn. P)dl 


repre- | Since 
:) Pr) 
(dl) dn dl . dn (6a) 
on d 
and 
d dO 4 
(dn) dl dn. al (6b) 
w yn 
| (see Fig. 3). 
| | PY d a d VY 0 | 
‘, dl . in +r ; ll.d 
this problem | e | yn (r yn ): ' yn” MM " 
) d~ 
= ldn. ? (¢ P)dl +7 Po” dn. dl} 
| M yn | 
,d/_d) Vd Lr »4\,. 
— (rP) +rl ) 
K,,,| au am 2 1} dn | (7 








Fig. 3. Enlarged diagram of the element a, b, c, d. 


CONTINUITY 
The rate of flow of fluid through the elementary 
ring V.2ardn. 
Since the rate of flow between two streamlines 


is constant, 
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\ 


hh 


.2ardn) 0 


dir) 


Nl 


‘ 


dn —rV 


From equation (6b), 


dirV) 


ow 


dn 


d(rV) 
MM 


Bounbary CONDITIONS 


The flow must be symmetrical with respect to 
adjacent particles and the velocity at the surface 
of each particle is zero 


RESULTANT SHEARING Force ON SPHERE 


Consider a ring of length dil at the surface of the 
sphere, as shown in Fig. 4. The shearing force 
on the ring can be resolved into a vertical and 
the net horizontal 


a horizontal component ; 


component must be zero. 


Fig. 4. Diagram of elementary ring on the surface of a 


particle. 


The vertical force on the ring 


dh 


rT<a 


=2radly 2Y) cos 0 = 27 an(>* | 
ON Sema on 


(since di cos @ = dh). 


.. resultant force on sphere 


b 
[29 ay 


. 


-b 


ot dh (9) 


on e 


In order to simplify the above equations and 
to enable a solution to be obtained, it will be 
assumed that d)P/dh is constant at a given 
horizontal plane at all points between r = a and 
r= R. This is equivalent to taking the value 
of the angle @ between the streamlines and the 
vertical as very small; this will not be appre- 
ciably in error, except near the surface. We 
therefore take 

cos 9 = 1 


Equation (7) then becomes 


d , av 
p 
on on 


, 292) 


wr ; 
| ay = sin a) 
. * > 
ot Jan = | r ol dn 
ATT) im M 
f(r, h) 
3P dP dh 
Mw MM 
~ oP 
~ Dh 


= ccs Ow Rend | = 
Ay 


. 
a 
. (r 
on 


Now since P 


dh 


(taking 7 sin @ which is 


small). 
Thus from equation (10), 


dV | 
r =-Iif 
on # d 
(since 


1P F 


= cos ~1) 


+k 





~ gdh’ 2 
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(where & is an integration constant). 


dV 1dPr ik 


: on me fe dh 2 ’ 
. av 
When r = R, 0 
on 
1)P R? 
k 
pp dh 2 
and 
av 1 \P ‘ed 
on Qin Wh r 
Since 
vr 
cos8 = 1 
on 
av 1 OP k? 
: (r (12) 
»” 2p wh ’ 


rte di Rin) k 


l MP (" 


» [p2 
~&’ - 1 wie R? In a) 
2p dh \2 
1 4 3 ae 


» 


Rin” ) (13) 
a 


; 2 dh 


The mean velocity V, of flow of tluid based on the 
total area of the cylinder of radius R is given by 


cos @.2ardr 


cr=R 
n KR? V, = | 


J r=a 


Substituting for V from equation (13) and 


putting cos @ = 1, 
7k VY, 
[’ ®¥ 1 OP (Pr —@ 
Jr=<e 2p dA 2 
rok 


9 aP 2 _ Pa a we (7, )(In = )| 
2 ¢ x 


2rrdr 


Rin’ 
é 


1 


pls ' 
MP pl, 
=? : = (14) 
yh Boel, 1/4 ¢ *(° = 4, R| 
R®) 1. ak In- 
8 | +3 rk) 3\rk) 3 "al 


The force on the sphere 


: on dV 
I = [0 eran(e ) 


. 


dh (equation 9). 


a 
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From equation (11), we have 


av 1 oP = 
An 2p dh (2 | a 


“| ' rR joe Mz) || dh 


5 lah kh 
+8 a\* 
7V. “| | Ur) )ah 
lja 4, a\* on R 
“. i P| | n) a 
; 2 
: (;,] th 
2 rVano){(. lja\* 4/a\? 4) R) > 
4 ‘| R| { a H : a 
. - » 
{38 R dhs 
taiabhadS [ lja\* 4/a\? 4 Ri 
J 14 | : | = * 
> alr) alr) 3a 
= 3a p Sf d(B,) (16) 


where £8, is a theoretical correction factor for use 
in conjunction with Stokes’ Law; its value is 
dependent only on the concentration. 


Table 1. y as a function of a, R 








a/R 

0 0 
0-001 O-1218 
0-005 0-1651 
0-01 0-1947 
0-035 O-2717 
O05 00-3328 
1 O-4752 
O-2 O-7987 
O38 1-2599 
O-4+ 1-9656 
0-5 31819 
0-6 54566 
O-7 10-6136 
0-8 25-597 
0-9 143-57 
0-95 832-7 

1 




























J. F. Ricnarpson and W. N. Zant: The sedimentation of a suspension of uniform spheres 


YVoueue—e 





© O05 O! O'S O2 O25 O3 O35 O4 








Fig. 5(a). 
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$05 O55 O68 


Fig. 5(b) 


y as a function of a/R (3 graphs.) 








In order to facilitate the graphical evaluation 
of the integral in equation (15), values of y are 
calculated for various values of a/R (Table 1) 
and plotted (Fig. 5). 


Here 


(5) - 


; | *) - : (“)- : In 
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Table 2 








- 
bR = 11267C 
configuration I 


b/R = 12254 C 
| configuration II 
| 






c | RB 








¢ 
ceded R(Il) B, b R (1) B, 

= | | 
0-05 1-271 | 0-274 0-78 O-415 1-32 
0-1 1-634 | 0-388 1188 = 0-523 211 
O15 2-128 0-474 1-68 0-599 2-99 
0-2 2-818 0-548 2-335 0-659 | 15 
0-25 3806 0-612 3-21 O71 | 5-62 
O-% 5248 «0-671 407 =O-754 7:58 
0-35 7-413 | 0-725 6-36 0-74 10-5 
O-4 10-73 | O-T775 8815 0-83 151 
0-45 15-98 | O821 138 0-864 | 2h4 
0-5 25-12 0866 24 37 O-804 38-5 
0-55 40-83 0-909 51°27 0-921 6 
0-606 75-86 00-9525 190 09525 190 


























O3 O4 O5 


~_ 





Fig. 6. 


Experimental or theoretical correction factors vs. 
concentration 


© Expermental correction factor vs. concentration, 

A Theoretical correction 
concentration. 

© Theoretical correction fector (configuration I) vs. con- 
centration, 

4+. Theoretical correction factor (BURGERS) vs. concentra- 
tion. 


factor (configuration II) vs. 
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In Table 2, values of 6 R and 8, are given for 
the two configurations as a function of the 
volumetric concentration C, and 8, is plotted 
against C in Fig. 6. 


A sample calculation is given below in Table 3. 
Evaluation of the theo- 


Table 3. Sample calculation. 

















retical correction factor for b/R = 0-548 
hob 1—(h/by @ Ji (*\* -- oe Y 
b b R Rb 

0 1 1 0-548 4:12 
Ol Oooo 0-994909 0-545 4-09 
Oz 0-06 0-9798 0-536 3-87 
O--4 0-84 0, 91652 0-502 3-23 
O-6 0-64 Os 0-439 2-37 
O-8 0-36 0-6 0-329 1-43 
oo 0-19 0-4558 0-239 0-96 
0-99 0-0199 0-14107 0-O774 0-41 
1 0 0 0 0 





co ES 4 
| | 














2 + 
y 

1} + +———— — 
O/R = 0-548 
AREA = 2-625 
Be = 8/9 * 2-625 =2-335 

l l | 
oO 0-2 O-4 0-6 o8 Oo 


hh=- 


S. 
b 


“1 
Fig. 7. Graphical evaluation of | y 
0 
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y is plotted against hb in Fig. 7, and the area 
under the curve is 2-625. Hence the theoretical 
correction factor = 8/9 x 2-625 = 2-335. 


EXPERIMENTAL VERIFICATION 


The results of a study of the sedimentation and 
fluidisation of particles of various shapes, in- 
cluding spheres, has been reported elsewhere 
[4]. For conditions of viscous flow, a relation 
was obtained, connecting the terminal falling 
velocity of a single spherical particle in an infinite 
medium (V,) with the velocity of fall (V,) of 
suspension of uniform spherical particles con- 
tained in a vessel where diameter was large 
compared with that of the particle. 


where « = the voidage of the suspension ( 


Vy (1 — cys (18) 
For a single particle settling at its terminal falling 
velocity, we have 

nd® 


F = 38ap Vid 
6 


(p, 


d? (p, 
Su 


pg 


and for a particle in a suspension 


F Oo V d B, - (p, 


pz 
alle, — v8 _ Ne (20) 
Is B Pe 8, 
where 8, is an experimentally determined cor- 
rection factor from Stokes’ Law. Comparing 
equations (18) and (20). 
B,=(1 —Cc)** (21) 
Values of 8,. are plotted against the volumetric 
concentration C, in Fig. 6, and are seen to lie 
between the two theoretical values. It would 
therefore appear that the configuration of par- 
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ticles is intermediate between the two theoretical 
arrangements and that arrangement 1 is 
approached more closely, especially at high 
concentrations. 

Also shown in Fig. 6, as a broken line, are the 
results of a theoretical analysis by Burcers [6], 
who investigated the effect of concentration on 
the settling velocity of a suspension of spheres. 
He concluded that for a random distribution of 
particles 
V, = V,(1 + 6875 Cy 


A (22) 

The agreement between Burcers’ theoretical 
expression and the experimental data is poor 
except at concentrations less than 0-2. The 
results of the present theory agree quite closely 
at high concentrations but extrapolate to give 
a correction factor of zero instead of unity at 
infinite dilution. This error appears to have a 
parallel in “ Stokes’ Paradox,” discussed by 
Birkuorr [5], in which the drag force exerted 
on a section of an infinitely long cylinder, with 
its axis parallel to the direction of flow, extra- 
polates to zero as the radius of the boundary walls 
approach infinity. 

The error in the present theoretical work 
appears to arise because an average pressure 
gradient )P/dh, has been used in the calculation 
of the drag force on particles of the suspension. 
At relatively high concentrations )P /dh is finite 
but at very low concentrations )P/dh tends to 
zero and it is not therefore permissible to sub- 
stitute for it. 

The theoretical treatment is, of course, not 
rigorous because no exact solution of the problem 
is possible with the mathematical techniques at 
present available, and it has been necessary to 
make assumptions about the distribution of the 
particles and the angles between the streamlines 
and the direction of flow. At high concentrations 
however a reasonably satisfactory correlation of 
the theoretical and practical figures has been 
obtained. 


ACKNOWLEDGMENT 
The authors are indebted to Professor D. M. 
Newitt, F.R.S., in’ whose laboratories the 


experimental work was carried out, for his 
constant help and encouragement. 


72 





J. F. Ricuarpson and W. N. Zaki: The sedimentation of a suspension of uniform spheres 


NOTATION radius of cylinder of fluid round a particle 
half distance between particles in hori- 
zontal direction 
distance from centre of sphere in a hori- 
zontal direction 
relative velocity between fluid and particle 
terminal falling velocity 
falling velocity of suspension 

= experimental and theoretical correction 
factors to Stokes’ Law 
density of the liquid 


- radius of circle at height (h) from centre 
of sphere = VY b? — h? 
= radius of sphere 
= volumetric concentration 
diameter of sphere 
= drag force 
= acceleration due to gravity 
= distance from centre of sphere in vertical 
direction = V b? — a? 


integration constants = density of spheres 
distance in the direction of streamlines — angle the streamlines make with the vertical 
= distance in the direction of equipressure 8 fl dh 
Y 
9/o0° b 


lines is given by the relation 8, 


pressure ML-'T-? porosity 
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Book reviews 


Silicones and their Uses. 
pp. 302 — XV. 


Ros Roy McGrecor : 
McGraw-Hill Book Co., New York, 1954. 
Price $6.00. 


This book is intended as a practical manual for engineers, 
metallurgists, chemists, etc. who intend to use silicones. 
There is little emphasis on the chemistry of the compounds, 
or on the reasons why silicones display such unusual 
properties, and non-chemical language is used wherever 
possible. 


An excellent introductory account of the early history 
of silicon chemistry leads, via a note on the commercial 
development of silicones, to the second, and most useful, 
of the five chapters of the book. Over 150 pages are 
spent in describing the various types of silicones that can 
be obtained — silicone fluids, resins, lubricants, rubbers, 
ete. The physical properties and applications of these 
materials have been obtained mainly from industrial 
pamphlets ete. issued by several American companies, 
particularly the Dow-Corning Corporation. The informa- 
tion, which is neatly summarised and often tabulated, 
constitutes a very valuable compendium. ‘The markedly 
non-critical approach is somewhat irksome, however, and 
there is a tendency to make silicones sound like a panacea 
for all ills. 


The section on silicone fluids covers such topics as 
stability in air, resistance to heat or chemicals, effect on 
organic rubbers, solubility in various solvents, effects with 
metals, boiling and freezing point, viscosity and shear 
stability, lubricating and dielectric water 
repellency, surface tension, specific heat, uses as anti- 
foaming agents, moulding, film, paper, textil> additives, 
and hydraulic fluids. 
type “continued exposure at high temperature ~ 
encountered, but on the whole the information provided 
under each heading is adequate and often excellent. 


properties, 


Occasional vague phrases of the 
are 


Dimethyl] silicones to which a few percent of a finely 
divided silica have been added have been termed * com- 
pounds,” and the properties and uses of such materials, 
considered in the next short section of chapter 2, include 
antifoaming agents, tap greases, and electrical applica- 
tions. This is followed by a surprisingly short section on 
silicone lubricants, and a particularly useful section on 
coating, laminating, release, water-repellant, moulding, 
and electrical resins. The following section on silicone 
rubber collects together for the first time information 
about chemical and physical properties (flexibility, hard- 
ness, tensile strength, elongation, compression set, 
chemical, oil, abrasion resistance, etc.), and includes a 
short account of the preparation and curing of the rubber. 


A short chapter on the physiological effects of silicon 
compounds is used mainly to emphasise the inertness of 
silicones ; it is followed by a chapter listing the uses of 
silicones in a variety of industries. The extent to which 





silicones are actually used as, say, * a defoamer for asphalt 
applied to highways ” as compared with, say, “ polishes 
for automobiles, furniture, and glass is not indicated, 
and many of the * uses ” listed would be more correctly 
classed as “ suggestions.” The last chapter, intended 
mainly for the non-chemist, gives a clear summary of the 
chemistry of silicon compounds and the methods used for 
the preparation of silicones. 


The bibliography of 152 references is disappointingly 
small when one considers the vast number of papers which 
appear annually on the subject. The book is well set out 
and printed, has few factual or typographical mistakes, 
and has a good index. It is essentially a non-crit’cal 
compilation of information and at times tends to read 
like a manufacturer's pamphlet. Nevertheless it fulfils 
a very useful purpose and is to quite an extent comple- 
mentary to Rocnow’'s Introduction to the Chemistry of the 
It can be recommended to those likely to use 
or consider silicones as a means of solving specialised 
R. N. Haszecpine. 


Silicones. 
problems. 


G. A. Morris and J. Jackson : Absorption Towers. 
Butterworths Scientific Publications, in association with 
Imperial Chemical Industries Limited, 1953. 159 pages. 


Price 30s. 


This book i; concerned primarily with the design of packed 
towers for the abs rption or desorption of gases. The 
emphasis is on the methods and data needed by an 
engineer to design full-scale towers which are technically 
and economically satisfactory. The methods described are 
based largely on the results of work carried out by Imperial 
Chemical Industries Limited, much of which has previously 
not been published. Since the aim of the book is different 
from that of a text-book, the fundamental theories of 
diffusion and analogies between fluid friction, heat trans- 
fer and mass transfer have been omitted, whereas con- 
siderable attention has been paid to economic factors and 
the design of liquid distribution systems. 


Both absorption coellicients and transfer units are used 
in various illustrative examples. Simple types of laboratory 
absorption equipment, such as the falling-film tower and 
the dise tower, are considered first. These towers are then 
used as standards with which the performance of packed 
towers can be compared. Gas-film coeflicients, for example, 
are estimated from data obtained in the falling-film tower, 
corrected by a gas-film packing factor. Liquid-film co- 
efficients are estimated from data obtained in the standard 
dise tower, corrected by a liquid-film packing factor. These 
individual coefiicients are combined in the usual way to 
give an overall coefficient which can be used for design. 
Individual and overall transfer units are estimated in a 
similar fashion. Systems in which absorption accompanied 
by chemical reaction occurs are discussed only very 
briefly. 










Several clearly worked examples are included to 
illustrate both the design of packed towers and the pre- 
diction of the performance of existing towers. There is a 
good discussion of wetting and flooding rates, pressure 
drop and economic gas rates for various packings. Several 
useful correlations and graphs are presented. Liquid and 
gas distribution systems are considered in detail. 


This book should be most helpful as a reference work 
for design engineers. It is well written and contains much 
useful information. Raymonp W. SourHworrtn, 


Inorganic Syntheses. Volume IV. Editor-in-chief, 
J. C. Battar, Jnr. McGraw-Hill Book Company, Inc., 
London. 1953. Price 36s. 


The aim of Inorganic Syntheses is to provide in book form 
complete and independently tested directions for the 
laboratory preparation of a large number of important 
inorganic substances, most of which are not obtainable 
commercially. Each volume contains a selection of 
compounds of elements distributed throughout the periodic 
table ; the present volume includes lithium nitride, silver 
(II) oxide, calcium hydrogen phosphates, iodochlorosilanes, 
nitrosyl and nitryl chlorides, diethyl phosphite, sodium 
superoxide, addition compounds of chromium (VI) oxide, 
potassium ferrate, and several anhydrous metal halides. 
Cumulative subject and formula indexes cover Volumes 
I-IV and greatly facilitate search for information. On 
the experimental side this volume fully maintains the high 
standard of its predecessors ; some of the descriptions of 
properties, however (e.g., a very elementary account of 
the classification of halides) might well have been omitted 
in a book written mainly for research workers in inorganic 
chemistry. A. SHARPE. 


C. K. Incotp : Structure and Mechanism in Organic 
Chemistry. G. Bell & Sons Ltd., London, 1953. vii + 
828 pp., 47 figures. Price £3. 17s. 6d. 

This work is highly imposing and offers much more than 
an ordinary text-book. It presents a picture of the 
development of theoretical organic chemistry, which may 
be considered historically very valuable in view of the 
large contribution made to this development by the 
author. Besides, the book contains a wealth of informa- 
tion which, no doubt, will prove a great help to countless 
research chemists in studying their problems. 


As a text-book for students of technology it is, 
undoubtedly, too exhaustive; the advanced student 
specializing in organic chemistry, however, will find it a 
storehouse of knowledge. 


In the initial chapters the author gives a clear picture 
of the development of the present structural formulation 
of organic compounds, Next he shows the reader how 
the modern theory is based on the physical study of 
matter and explains the terminology used in the following 
chapters. 
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In these chapters, as in the preceding ones, the matter 
discussed is amply illustrated with experimental data. 
Justly, a separate chapter has been devoted to the 
aromatic nucleus and the concept of the aromatic character. 
After this the author proceeds to give a classification of 
reagents and reactions. 

The nucleophilic substitution is extensively discussed, 
which is not to be wondered at, in view of the large 
contribution to this field which has been made by Prof. 
IncoLp and his school. 

In the last few years many books treating of theoretical 
organic chemistry have been published. In our opinion, 
however, the appearance of a work such as this, which, in 
its elaborateness, proves to be the fruit of many years of 
development, is a rare occurrence. 

With gratitude we accept this book, which enables us 
to survey with the author, the achievements in this field 
of organic chemistry. Dr. G. pe Jone, 


Rotanp 8S. Younc: Industrial Inorganic Analysis. 
London, Chapman and Hall Ltd., 1953. 
12 figs., 7 tables. Vrice 36s, 


viii + 368 pages, 
It has been the intention of the author — in view of, for 
example, the A.S.T.M. Methods of Chemical Analyses of 
Metals — to give a description of such methods of analysis 
as may occur in inorganic laboratory practice and which 
have not yet been described in full in any other book. It 
is a real pleasure to note that at long last there has 
appeared a book which, unlike practically all other 
analytical handbooks, does not start with 50-100 pages 
on general subjects which may well be considered as 
common knowledge. 

Forty-four elements arranged in alphabetical order are 
discussed. For the determination and separation of each 
element, the author presents several methods supple- 
mented with — on the whole — sound comments. 

Whereas the methods given for the analysis of some 
elements (Cr, Ni) are very good, it must be remarked that 
several others seem to be rather obsolete. For example, 
the book contains quite a large number of determinations 
in which use is made of hydrogen sulphide, although, as 
is well known, the presence of H,S in these tests is anything 
but desirable. 

The chapter on Miscellaneous Analyses and Data, 
contains about 35 pages dealing with the determination 
of acids in some special samples, gas analyses, examination 
of water, examination of xanthates, standard acids and 
bases and indicators. Partly because of its conciseness, 
this appendix is much less successful than the part on 
the analysis of the elements. 

The reference literature cited is chiefly by American 
authors. Much care has been given to the author and 
subject indexes. 

Apart from a few inaccuracies, the book provides much 
useful information to the more experienced analytical 
chemist. H, W. Dernum, 
































Letter to the Editors 


Layer formation in a horizontal extractor 


In an attempt to study mass transfer between two 
horizontal flowing viscous streams of carbontetrachloride 
with acetic acid as a solute and water, we met with the 
remarkable effect that no steady state conditions could 
be established even over very prolonged periods. The 
dimensions of the extractor used were 100 x 5 «x 5em 
(length, width and height), the streams entered the 
apparatus without undue disturbance and ran counter- 
currently. At a volume rate of 3 l/hour for each stream, 
even after 16 hours, therefore after a passege of 38 times 
the volume of the apparatus, it appeared that the acetic 
acid content of the carbontetrachloride had decreased 
20-30%, more than the increase in the acetic acid content 
of the water. The material balances were taken from 
the entering and leaving streams. 


It turned out that this very slow establishment of 
steady state conditions is caused by a strong tendency 
for layer formation caused by density differences of 
solutions of different acetic acid content. 


Proof is given by replacing during an extraction experi- 
ment the solutions of the otherwise same two liquids by 
coloured solutions of the same composition. It appeared 
that the colourless (upper) water layer was slowly replaced 
by the coloured water layer with the exception of a small 
zone just above the interface carbontetrachloride-water. 
For a volume rate of 3 l/hour of both streams at the 
uncoloured zone was 
A very 


entrance side of the water this 
4mm thick, at the exit side 7 mm, for hours. 
sharp separation between the two water lavers existed. 





Apparently, just above the interface a more or less 
stagnant layer of higher density is formed, which accumu- 
lates part of the acetic acid passing the interface. This 
indeed could be shown with a universal indicator. In- 
creasing the volume rate of the water stream was accom- 
panied by a decrease in the thickness of the layer as 
could be expected. If a stream of a coloured 5%, brine 
solution of density 1-04 g/cm" was introduced after the 
experiment described above the very spectacular observa- 
tion was made that the brine stream entered the colourless 
water layer like a wedge and split this colourless layer 
into two colourless layers. The layer near the interface 
of higher density than 1-04g/em® had a thickness of 
1-5 mm. 


For the carbontetrachloride stream the effect was the 
opposite, the coloured stream flowed along the interface. 
By a decrease in the acetic acid content, the density 
increases and in the direction of flow a mixing effect in 
a downward vertical direction was observed. 


It is obvious that the effects described strongly will 
hinder mass transfer measurements in a_ horizontal 
apparatus. Although the influence of the concentration 
on the density is well known, the effect on the behaviour 
of a horizontal laminar flowing liquid is greater than we 
expected. Because this also might be the case with 
colleagues working in the same field, we thought it worth 
while to bring our observations to your notice. 


C, Hapicutr 
P. M. Heertses 
H. TatsmMa 


University of Technology, 
Laboratory for Chemical Engineering, 
Delft, Holland. 





